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ABSTRACT 

In the present paper, we study how the dark matter density profiles of dwarfs galax- 
ies in the mass range 10 s — 10 10 M Q are modified by the interaction of the dwarf in 
study with the neighboring structures, and by changing baryon fraction in dwarfs. 
To this aim, we determine the density profiles of the quoted dwarfs by means of Del 
Popolo (2009) which takes into account the effect of tidal interaction with neighbor- 
ing structures, those of ordered and random angular momentum, dynamical friction, 
response of dark matter halos to condensation of baryons, and effects produced by 
baryons presence. As already shown in Del Popolo (2009), the slope of density profile 
of inner halos flattens with decreasing halo mass and the profile is well approximated 
by a Burkert profile. We then treat angular momentum generated by tidal torques 
and baryon fraction as a parameter in order to understand how the last influences the 
density profiles. The analysis shows that dwarfs who suffered a smaller tidal torquing 
(consequently having smaller angular momentum) are characterized by steeper profiles 
with respect to dwarfs subject to higher torque, and similarly dwarfs having a smaller 
baryons fraction have also steeper profiles than those having a larger baryon fraction. 
In the case tidal torquing is shut down and baryons are not present, the density pro- 
file is very well approximated by an Einasto profile, similarly to dwarfs obtained in 
dissipationless N-body simulations. We then apply the result of the previous analysis 
to the dark matter halo rotation curves of three different dwarfs, namely NGC 2976, 
known to have a flat inner core, NGC 5949 having a profile intermediate between a 
cored and a cuspy one, and NGC 5963 having a cuspy profile. After calculating baryon 
fraction, which is ~ 0.1 for the three galaxies, we fitted the rotation curves changing 
the value of angular momentum. NGC 2976, has an higher value of ordered angular 
momentum (A ~ 0.04) with respect to NGC 5949 (A ~ 0.025) and in the case of NGC 
5963 the very steep profile can be obtained with a low value of A (A ~ 0.02) and 
also decreasing the value of the random angular momentum. In the case of NGC 2976 
tidal interaction with M81 could have also influenced the inner part of the density 
profile. We finally show how the inner density profile correlates with Karachcntsev et 
al. (2004), and Karachcntsev and Kashibadze (2006) tidal index for dwarfs and LSBs. 

Key words: cosmology-theory-large scale structure of Universe-galaxies-formation 



1 INTRODUCTION 

Notwithstanding dwarf galaxies are the most diffused type 
of galaxies in the Universe, they have received less atten- 
tion than spiral and elliptical galaxies. Their relevance is 
noteworthy in cosmology since their study may bring new 
insights in the processes present in the early universe and 
since the dwarfs that we observe today are probably the sur- 
vivors of the building blocks of larger galaxies. Dwarfs can be 
classified in dwarf spiral galaxies, blue compact dwarf galax- 
ies (BCDs), dwarf irregular galaxies (dlrrs), dwarf elliptical 

* E-mail: antonino.dclpopolo@unibg.it 



galaxies (dEs), dwarf spheroidal galaxies (dSphs), and tidal 
dwarf galaxies. They are all characterized by small size, low 
mass, low metallicity, and low luminosity, and are moreover 
characterized by an old population (Grebel 2001). The so 
called morphological segregation is the difference observed 
in spatial distribution of dSphs, DEs, and dlrrs, and is prob- 
ably originated by the effects that environment has in the 
evolution of galaxies. As previous reported, dwarfs are im- 
portant cosmological probes and it is not an incident that 
the missing satellites problem (Klypin et al. 1999; Moore et 
al. 1999) and the cusp-core problem, (de Blok 2010), namely 
the discrepancy of observed and simulated density profiles 
in the inner kpc of the halo, are to them connected. For 
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what concerns the second problem, on which we shall con- 
centrate in the rest of the paper, comparison of rotation 
theoretical curves, obtained by Moore (1994) and Flores & 
Primack (1994), with HI rotation curves of late-type, gas- 
rich, dwarf galaxies lead to the conclusion that these objects 
are dark matter dominated and their halos are best fitted 
by constant-density core, instead of the theoretical predicted 
cuspy profiles. This lead to the conclusion that CDM pre- 
dictions on small scales are incompatible with observations 
(Carignan & Beaulieu 1989; Broeils 1990) of dwarf galaxies 
(Moore 1994). Also LSB galaxies seem to indicate that the 
shape of the density profile at small scales is significantly 
shallower than what is found in numerical simulations (e.g., 
de Blok, Bosma, and McGaugh 2003), and a similar indica- 
tion is given by spheroidal galaxies (e.g., Ursa Minor (Kleyna 
et al. 2003), Draco dSph (Magorrian 2003)). 

According to the studies of the density profile of dwarfs, 
the inner part of density profiles is characterized by a core- 
like structure (e.g., Flores & Primak 1994; Moore 1994; 
Kravtsov et al. 1998; de Blok & Bosma 2002 (hereafter 
dBB02); de Blok, Bosma & McGaugh 2003; Gentile et al. 
2004, 2006; Blaise-Ouelette et al. 2004, Spano et al. 2008, 
Kuzio de Naray et al. 2008, 2009, and Oh et al. 2010 (here- 
after OhlO)), but some studies found that density profiles 
are cuspy. For example, Hayashi et al. (2004) found that 
about 70% of galaxies in the sample they used have rotation 
curves consistent with steep slopes. Other studies concluded 
that density profiles are compatible with both cuspy and 
cored profiles (van den Bosch et al. 2000; Swaters et al. 2003; 
Spekkens et al. 2005; Simon et al. 2005 (hereafter S05); de 
Blok et al. 2008 (hereafter dB08)). In the case of NGC 2976, 
4605, 5949, 5693, 6689, S05 found that the power-law index, 
a, describing the central density profile, span the range from 
a = (NGC 2976) to a = 1.28 (NGC 5963), with a mean 
value a = 0.73, shallower than that predicted by the sim- 
ulations, and scatter in a from galaxy to galaxy is equal 
to 0.44, which is 3 times as large as in cold dark matter 
(CDM) simulations. dB08, by using the THINGS sample, 
found that for galaxies having Mb < — 19 a NFW profile or 
a pseudo-isothermal (PI) profile statistically fit equally well 
the density profiles of galaxies, while for Mb > —19 the 
core dominated PI model fits significantly better than the 
NFW model. In other terms, for low mass galaxies a core 
dominated halo is clearly preferred over a cusp-like halo, 
while for massive, disk dominated galaxies, the two models 
fit apparently equally well. 

The previous arguments lead to some important ques- 
tions: is there a correlation between the value of slope of 
galaxies tidally influenced by external galaxies and the tidal 
field? How much the baryons in dwarfs influence their fi- 
nal density profiles? Since some studies (e.g., S05, dB08) 
hint to a possible correlation between inner slope of dwarfs 
and their mass and interaction with neighbors, would be 
important to verify from a theoretical point of view if the 
quoted correlations are admissible. On the other hand, the 
morphology-density relation illustrates that environment is 
an important factor in dwarfs evolution. 

This paper has a double goal: the first one is to see 
how changes in baryons fraction of galaxies and magnitude 
of tidal interaction modify the density profile of galaxy. It 
is important to stress that in Del Popolo (2009) the effects 
of baryon presence was already taken into account, but we 



used fixed values of baryon fraction, according to the type 
of galaxies studied, and angular momentum was obtained 
and fixed to the value given by the tidal torque theory 
(TTT). While nowadays we are able to compute with rea- 
sonable precision the evolution of dark matter component, 
physical processes driving galaxy formation, like the evolu- 
tion of baryonic matter, cooling, star formation, black hole 
formation, and SN and AGN feedback, are extremely com- 
plex and still debated. However, these processes, usually re- 
ferred to as sub-grid physics, can be taken into account with 
new computer models. One example is the hydrodynamical 
simulation presented by Governato et al. (2010) (hereafter 
G10), which obtained bulgeless dwarf galaxies with shallow 
inner dark mater profile. In their study, the authors took 
account of baryons and showed how supernovae explosions 
affect galaxy evolution. The result obtained needed a inten- 
sive use of supercomputers for long CPU times. Analytical 
and semi-analytical models, like the one presented in Del 
Popolo (2009), have several advantages when compared to 
simulations: a) flexibility (one can study the effects of phys- 
ical processes one at a time including subgrid processes); b) 
computational efficiency (it takes about 10 s to compute the 
density profile of a given object at a given epoch on a desk- 
top PC, Ascasibar et al. 2007). If we compare our model 
to G10 simulation, as we will show in the following, we ob- 
tain accurate density profiles, but at the same time we can 
obtain different density profiles for different values of angu- 
lar momentum and baryon fraction, with small CPU times. 
Studying how baryon fraction or tidal interaction with large 
neighbors influence the density profile, by means of G10 sim- 
ulations, would imply to perform each time a new simula- 
tion, with a prohibitive computational cos10. 

The second goal of the present paper is that of using 
the results of the first goal (changes of density profiles with 
baryon fraction and angular momentum) in order to clarify 
if the density profile of some observed galaxies with cuspy, 
cored, and intermediate density profiles can be re-obtained 
in our model. So, in the present paper, we analyze how dif- 
ferent tidal torques and baryons content of haloes influence 
dwarfs density profile. 

The paper is organized as follows: in Section 2, we sum- 
marize the model used. In Section 3.1, we determine the 
dwarfs density profiles when tidal torques and baryons frac- 
tion is changed. In Section 3.2, we compare the previous re- 
sults with three different galaxies, namely NGC 2976, NGC 
5949, and NGC 5963, we discuss why the three dwarfs have 
different profiles, and finally extend the study to some other 
galaxies. Section 4 is devoted to discussions and conclusions. 



2 SUMMARY OF DEL POPOLO 2009 MODEL 

In the present study, the dwarfs galaxies dark matter haloes 
are formed using the analytical method introduced in Del 
Popolo (2009), which is summarized in the following. We 
added a discussion on the way baryon fraction was calcu- 
lated, not present in Del Popolo (2009) (hereafter DP09). 
The spherical infall model (SIM), was studied for the 

1 For this reason, in Governato et al. (2010) simulation were fun- 
damentally calculated two isolated haloes with a fixed value of 
baryon fraction. 
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first time by Gunn & Gott (1972) and widely discussed and 
improved in a large series of following papers (Fillmore & 
Goldreich 1984 (hereafter FG84); Bertschinger 1985; Hoff- 
man & Shaham 1985 (hereafter HS); Ryden & Gunn 1987 
(hereafter RG87); White &Zaritsky 1992; Avila-Reese et al. 
1998; Nusser & Sheth 1999; Henriksen & Widrow 1999; Sub- 
ramanian et al. 2000; Del Popolo et al. 2000; Hiotelis 2002; 
Le Delliou & Henriksen 2003; Ascasibar et al. 2004; Williams 
et al. 2004). In SIM one considers a spherical overdensity in 
the expanding universe, that behaves exactly as a closed 
sub-universe because of Birkhoff 's theorem, and studies its 
evolution. A bound mass shell of the quoted overdensity will 
expand from an initial radius Xi to a maximum one, x m , 
which is usually termed apapsis or radius of turn-around, 
x ta (Peebles 1980): 



= g{xi) = Xi 



1 + S, 



(1) 



where Qi is the density parameter at time ti, and the mean 
fractional density excess inside the shell, as measured at cur- 
rent epoch to, assuming linear growth, can be calculated as: 



5i 



S(y)y dy 
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Eq. ((2)1 furnish a relation among the final time averaged 
radius of a given Lagrangian shell and its initial radius. 

The initial overdensity 5i(xi) can be calculated once a 
spectrum of perturbations is known. Throughout Del Popolo 
(2009), we adopted a ACDM cosmology with WMAP3 pa- 
rameters, n m = 1 - f2 A = 0.24, Q A = 0.76, Sl b = 0.043 
and h — 0.73, where h is the Hubble constant in units of 
100 km Mpc^ 1 . The variance as = 0.76 (Romano-Diaz 
et al. 2008) of the density field convolved with the top hat 
window of radius 8 A -1 Mpc -1 was used to normalize the 
power spectrum. 

Expressing the scaling of the final radius, x, with the 
initial one by relating x to the turn around radius, x m , the 
final radius can be expressed as: 

x = f(xi)x m (3) 

being / the collapse parameter. Using mass conservation, 
and assuming that each shell is kept at its turn-around ra- 
dius, the shape of the density profile is given by (Peebles 
1980; HS; White & Zaritsky 1992): 



Pta{Xm)= Pi{Xi){—) -T— 



(4) 



After reaching maximum radius, a shell collapses and 
will start oscillating and it will contribute to the inner shells 
with the result that energy is not an integral of motion any- 
more and the collapse factor, /, is no longer constant. The 
following evolution of the system can be described assum- 
ing that the potential well near the center varies adiabat- 
ically (Gunn 1977, FG84). If a shell has an apapsis radius 
(i.e., apocenter) x m and initial radius Xi, then the mass in- 
side Xm is obtained summing the mass contained in shells 
with apapsis smaller than x m (permanent component, m p ) 
and the contribution of the outer shells passing momentar- 
ily through the shell x m (additional mass m a dd), and can be 
expressed as: 

4 3 - 

m T {xm) = m p (xm) + m add (x m ) = -np bA x l (1 + Si) (5) 
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where p b> i is the initial background density, R is the radius 
of the system (the apapsis of the outer shell) and the dis- 
tribution of mass m(x) — m(x m ) is given by Eq. @, while 
Px m (x) is the probability to find the shell with apapsis x in- 
side radius x m , calculated as the ratio of the time the outer 
shell (with apapsis x) spends inside radius x m to its period. 
This last quantity can be expressed as 
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where x p is the pericenter of the shell with apsis x and v x (rf) 
is the radial velocity of the shell with apapsis x as it passes 
from radius rj. 

This radial velocity can be obtained by integrating the 
equation of motion of the shell: 



dv r _ h 2 (r,u) + j 2 (r,u) 
dt ~ r 3 



„, . dr A 
GW ^+3 r 



(8) 



where h(r, v) is the ordered specific angular momentum 
generated by tidal torques, j(r,u) the random angular mo- 
mentum (see RG87 and Del Popolo 2009), G(r) the gravita- 
tional acceleration, A the cosmological constant and p the 
coefficient of dynamical friction. The final density profile is 
given by (Eq. A18, DP 2009): 



p(x) = 



pta{x m ) 



1 + 



din/ 



(9) 



As seen by the previous discussion, the previous model 
(i.e., Del Popolo 2009) takes into account several effects, 
such as ordered angular momentum, h, and random angular 
momentum, j, dynamical friction and adiabatic contraction 
of dark matter (AC). 

Ordered angular momentum is due to the tidal torques, 
r, exerted by a field of random density fluctuations on proto- 
halos (Hoyle 1953; Peebles 1969; White 1984; Ryden 1988; 
Eisenstein & Loeb 1995; Catelan & Theuns 1996), it can be 
calculated when the r.m.s. torque, r(r), is known, by inte- 
grating torque over time (see Del Popolo 2009 Appendix C, 
for details). Before going on, we recall that common practice 
is that of expressing the total angular momentum in terms 
of the dimensionless spin parameter 



A 



L\E\ 



i/-> 



GM^ < (10) 
where L is the angular momentum, summed over shells, 
and E is the binding energy of the halo. As several authors 
showed (e.g., Vivitska et al. 2002 ), the A parameter is well 
described by a log-normal distribution (e.g. Vivitska et a. 
2002), with a maximum of A = 0.035 and 90% probability 
that A is in the range 0.02-0.1 (Vivitska et al. 2002). While 
Catelan & Theuns 1996 showed a dependency of A on the 
peak heigh10 Knebe & Power (2008) and Cervantes-Sodi at 



<5(0)/cr, where a is the mass variance filtered on a scale 



R 



3 v is the central peak height given by v = <5(0)/cro(-R), where 8(r) 
is the overdensity in the radius r, and uq(R) the mass variance 
filtered on scale R.. 
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al. (2008a, b) (see Fig. 5 of Cervantes-Sodi et al. 2008a; and 
Fig. 3 of Cervantes-Sodi et al. 2008) showed a dependence 
on environment, cosmology, and mass, and Cervantes-Sodi 
et al. (2010b) found a correlation between A and cluster en- 
vironment. So, while the baryons in bright disc galaxies seem 
to have spin parameters similar to those of their host haloes, 
dwarfs and LSB disc galaxies may tend to be associated 
with a higher spin parameter. For example, van den Bosch, 
Burkert & Swaters (2001, hereafter BBS) studied the spin 
in a sample of 14 low-mass disc galaxies with an estimated 
average of virial velocity V ~ 60 km/s. They found an aver- 
age spin parameter about 50% larger than that of the dark 
haloes (see also Mailer & Dekel 2002, Fig. 8 and Dekel & 
Woo 2003). Also Boissier et al. (2003), showed in their study 
that 35% of all galaxies (in number) with 0.06 < A < 0.21 
are LSBs (see also Jimenez et al. 1998). According to several 
analytical papers (Nusser 2001; Hiotelis 2002; Le Delliou & 
Henriksen 2003; Ascasibar et al. 2004; Williams et al. 2004), 
objects with larger angular momentum should have flatter 
profiles, one therefore expects that these high spin objects 
should typically have shallower density cusps. Moreover the 
central parts of these galaxies are also much less dense than 
the simulations indicate. 

Random angular momentum (RG87, see Appendix C2 
of Del Popolo 2009), is connected to random velocities (see 
RG87), and in all previously quoted papers only this kind of 
angular momentum was taken into consideration. The usual 
approach (Nusser 2001; Hiotelis 2002; Ascasibar et al. 2003) 
consists in assigning a specific angular momentum at turn 
around: 

j oc VGMr m (11) 

With this prescription, the orbital eccentricity e is the same 
for all particles in the halo (Nusser 2001). We assigned ran- 
dom angular momentum to protostructures expressing the 
specific angular momentum j through the ratio e = ( ) , 
and left this quantity as a free parameter (Avila- Reese et al. 
1998) . This procedure is justified by N-body simulations of 
halo collapse: for CDM halos, N-body simulations produce 
constant < r r min >~ 0.2 ratios of dark matter particles in 
virialized haloes. Ascasibar et al. (2007), showed that there 
is a dependence on the quoted ratio on the dynamical state 
of the system: major mergers are well described by constant 
eccentricity up to the virial radius, relaxed systems are more 
consistent with a power-law profile. Minor mergers are in 
the middle. In the present paper, we use Avila-Reese et al. 
(1998) method with the correction of Ascasibar et al. (2003) 
(see their section 3.2 and their Eq. 32). 

Dynamical friction was taken into account by introduc- 
ing the dynamical friction force (see Appendix D of Del 
Popolo 2009) in the equation of motions previously de- 
scribed. 

We also took into account adiabatic contraction (AC) of 
dark matter halos in response to the condensation of baryons 
in their centers, leading to a steepening of the dark matter 
density slope. Blumenthal et al. (1986) described an approx- 
imate analytical model to calculate the effects of AC. More 
recently Gnedin et al. (2004) proposed a simple modification 
of the Blumenthal model, which describes numerical results 
more accurately. For systems in which angular momentum is 
exchanged between baryons and dark matter (e.g., through 
dynamical friction), Klypin et al. (2002) introduced a mod- 



ification to Blumenthal's model. The adiabatic contraction 
was taken into account by means of Gnedin et al. (2004) 
model, and Klypin et al. (2002) model was used to take also 
account of exchange of angular momentum between baryons 
and dark matter. 

For what concerns how baryon fraction and baryon dis- 
tribution was chosen, we want to recall that the model of 
Del Popolo (2009) contained both DM and baryons. In that 
paper and in its Appendix E, we discussed the quantity of 
baryons introduced in the model and how their initial mass 
distribution and the final distribution was fixed. Here, we 
speak again of baryons and baryon fraction, since here we 
study how the density profiles changes in terms of change of 
the baryon fraction. Till now, almost all N-body simulations 
on density profiles have dealt just with dark matter halos. 
However, we know that halos contain baryons, that influ- 
ence the evolution of dark matter. The effect of baryons on 
dark matter haloes was studied in El-Zant et al. (2001) and 
Romano-Diaz et al. (2008) simulation. Full SPH simulations 
have been performed recently by G10, showing the strong 
effect of baryon physics on the final dark matter profile. 

In the present paper, baryon fraction was calculated as 
follows. 

As discussed in the introduction, dwarf galaxies can be 
classified in galaxies having an high content of gas (e.g., 
dwarf irregular galaxies) and in galaxies having a low content 
of gas (e.g., dwarf elliptical galaxies or dwarf spheroidal). 
Moreover, baryon content changes in the galaxy from the 
core to the outskirts (see Hoeft et al. 2006, Fig. 5). The uni- 
versal baryon fraction, /;,, has been inferred by observations 
involving different physical processes (Turner 2002). Perci- 
val et al. (2001), by using the Two Degree Field Galaxy 
Redshift Survey, obtained a value of 0.15 ± 0.07 , Jaffe et 
al. (2001) obtained f b = 0.186±o:oJg 

combining several data 
sets, while by using WMAP collaboration data (Spergel et 
al. 2003) (Q b = 0.047 ± 0.006; On = 0.29 ± 0.07), one ob- 
tains fb — 0.16. However, the content of baryons in galaxies 
is usually much smaller than f b , since after baryons cool and 
form stars, the baryon-to-total mass at the virial radius is no 
longer equal the previous quoted universal baryon fraction 
and may deviate from it depending on feedback effects, hi- 
erarchical formation details, and heating by the extragalac- 
tic UVB flux. One has to take into account that for large 
galaxies, only about half of the baryons, in principle avail- 
able within the virial radius, goes into the central galaxy. For 
dwarf galaxies, this amount may be much smaller because of 
feed-back effects, hierarchical formation details, and heating 
by the extragalactic UVB flux. Geha et al. (2006) studied 
the baryon content of low mass dwarf galaxies and Hoeft et 
al. (2006) studied structure formation in cosmological void 
regions investigating to which extent the cosmological UV- 
background reduces the baryon content of dwarf galaxies. 
Hoeft & Gottlober (2010) studied also the case of cosmo- 
logical environment. As shown from McGaugh et al. (2010) 
(see their Fig. 2), baryon fraction deviates from the cosmic 
baryon fraction: while on dwarf scale only 1% of the ex- 
pected baryons are detected, on cluster scale most of them 
are detected. The detected baryon fraction is the ratio of the 
baryon fraction, F b = Mb/Msoo, and the universal baryon 
fraction, /;,: 

fd = (M b /M 500 )/f b = F b /f b (12) 
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where M500 is the mass contained within a radius where 
the enclosed density is 500 times the critical density of the 
Universe, and the baryon mass Mj = M* + M gas , where 
M» is the stellar mass and M gaa is that of gas. For example 
for a dwarf of M b = 2 x lO 7 Af , with M500 = 2 x 1O 9 M 
(see Table 2 in McGaugh et al. 2010), f d is ~ 0.04, and the 
baryon fraction Ft = Mb/M^oo = 0.01. 

The result of our model (extensively described in DP09) 
shows that the structure of the inner density profiles is 
strongly influenced by baryon physics and exchange of angu- 
lar momentum between baryons and dark matter. So, sum- 
marizing, in our analytical model the density profiles are: 

a) core-like, in the inner part of the density profile, when 
we take into account the presence of baryons. In the outer 
parts the result is the same of simulations. 

b) The results of dissipationless simulations, namely cuspy 
profiles, are re-obtained, with very good accuracy, eliminat- 
ing the baryons, not taken into account in dissipationless 
simulations. The results are in agreement with other mod- 
els, e.g. Mashchenko et al (2005, 2006, 2007), Romano-Diaz 
et al. (2008), El-Zant et al. (2001), and G10. 

As reported in the introduction and also in this section, 
the model of the present paper has some peculiar features 
that makes it unique: a) the other semi-analytical models 
take into account only random angular momentum, while 
the present paper takes also into account ordered angular 
momentum. Other studies, concentrate only on the study of 
one effect: angular momentum or dynamical friction or adi- 
abatic contraction, while we study their joint effect, b) We 
can take into account sub-grid physics accurately and with- 
out using too high CPU time, and this allow us to study how 
change of angular momentum or baryon fraction influence 
density profiles. 



3 RESULTS AND DISCUSSION 

3.1 Density profiles, torque and baryons fraction 

While very successful on large scale, ACDM model suffers 
from a number of problems on galactic scales (Flores & Pri- 
mack 1994; Moore 1994; Ostriker & Steinhardt 2003). One 
of these, the cusp-core problem is particularly troublesome 
and despite many efforts it has not entirely been resolved 
to date. In DP09, we showed how the flattening of the in- 
ner slopes of halos is produced by the role of angular mo- 
mentum, dynamical friction, and the interplay between DM 
and baryon component. In peculiar, we studied the density 
profiles of haloes in the mass range 10 8 — 1O 15 M . As previ- 
ously reported, in that study we moreover took into account 
the ordered angular momentum, due to tidal torques (see 
Section 2 of the present paper and Section C of DP09), ex- 
perienced by proto-halos in the neighborhood of the studied 
object. It would be then interesting to study the effects of 
different angular momenta and baryon fraction on the den- 
sity profile of a given halo. So, in the following, using DP09, 
we will generate density profiles of haloes in the mass range 
10 s - 10 10 A/ o and we will study the effect of different or- 
dered angular momentum and baryon fraction on the dwarf 
galaxy density profile. In a future paper, we also aim to 
study observationally the inner slopes of dwarfs who has 
no interacting companions, and those more close to parent 



galaxies. In Sect. 3.1, we use the baryons fraction of Mc- 
Gaugh et al. (2010). In Sect. 3.2, where we compare the 
results with three observed galaxies, we will discuss how the 
baryon fraction is chosen, for those peculiar cases. The result 
of the calculation are shown in Fig. la, in which we plot- 
ted the density profiles of dwarfs with mass 1O 8 M (dashed 
line), 1O 9 M (dotted line), and 10 10 A/ Q (solid line). This 
panel shows the results for haloes density profiles obtained 
in DP09 (which we will use as reference haloes). The angu- 
lar momentum in Fig. la was obtained by the tidal torque 
theory as in DP09 and, in the case of the halo of 1O 9 M has 
a value similar to UGC 7399 (see BBS), namely /i* ~ 400 
kpc km/s (A ~ 0.04), where h is the specific ordered angular 
momentum. The baryon fraction (fd, — 0.04)0 wa s fixed as 
described Sect. 2. The quoted plot shows flat density pro- 
files well fitted by a Burkert profile, whose functional form 
is characterized by: 

PW = (1 + r/r )[l + (r/r ) 2 ] (13) 

where p ~ p s and r ~ r a (El-Zant et al. 2001). 

A slight steepening of the density profile is visible in the 
case 10 10 M in which case the inner slope of the profile is 
a ~ 0.2, while a ~ for 1O 8 M -1O 9 M . 

As previous discussed in Sect. 2, the ordered angular 
momentum acquired by proto-structures depends on mass 
and is higher for dwarf galaxies, moreover the values of A 
has a range ~ 0.02-0.2 (BBS; Boisier et al. 2003). So, after 
calculating the density profiles for value given by the tidal 
torque theory (TTT), we considered the effects of changing 
the values of ordered angular momentum in the range of the 
allowed values. In Fig. lb, we repeated the same calculation 
of panel (a) changing the amplitude of the ordered angular 
momentum (which is proportional to the tidal torque) to 
2h, , leaving unchanged the baryon fraction to fd, . In this 
case, the profiles are, as expected, flatter with respect to the 
reference case, and for the same mass as the reference case 
(1O 8 M , 1O 9 M , 1O 1O M ) the slope is a = 0, in all three 
cases. In Fig. lc, we plot the density profiles for haloes hav- 
ing a value of angular momentum /i*/2, (1/2 that of the ref- 
erence case) , leaving unchanged the value of baryon fraction 
to fd, ■ In this case, the inner slope for haloes of the same 
mass (1O 8 M , 1O 9 M , 1O 1O M ) are, respectively, given by 
a ~ 0.4, a ~ 0.5, and a ~ 0.6. The previous results show a 
steepening of the slope with decreasing value of angular mo- 
mentum, in agreement with several previous results (Sikivic 
et al. 1997; Avila-Reese et al. 1998; Nusser 2001; Hiotelis 
2002; Le Delliou & Henriksen 2003; Ascasibar et al. 2004; 
Williams et al. 2004). The three panels show that: a) the less 
massive is a halo the less concentrated it is; b) larger mass, 
and higher angular momentum haloes are characterized by 
smaller inner slopes. 

Point (a) can be explained as follows: higher peaks (larger 
v), which are progenitors of more massive haloes (Peacock & 
Heavens 1990; Del Popolo k, Gambera 1996, Gao & White 
2007) have greater density contrast at their center, and so 
shells do not expand far before beginning to collapse. This 
reduces the angular momentum acquired and allows haloes 
to become more concentrated. Point (b) can be explained 

4 Note that the "*" in fd, and in h, does not "stand" for stellar, 
as in the case of M, 
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Figure 1. Shape changes of dark matter haloes with angular momentum. In panels (a)-(c), the dashed line, the dotted line and the 
solid line represent the density profile for a halo of 10 8 Mq, 10 9 Mq, and 1O 1O M0, respectively. In case (a), that is our reference case, 
the specific angular momentum was obtained using the tidal torque theory as described in DP09. The specific angular momentum for 
the halo of mass 10 9 Mq is /i* ~ 400 kpc km/s (A ~ 0.04) and the baryon fraction /j ~ 0.04. In panel (b) we increased the value 
of specific ordered angular momentum, h*, to 2h t leaving unmodified the baryon fraction to and in panel (c) the specific ordered 
angular momentum is /i*/2 and the baryon fraction equal to the previous cases, namely f^. Panel (d) shows the density profile of a 
halo of 10 Mq with zero ordered angular momentum and no baryons (solid line), while the dashed line is the Einasto profile. 



as folllows. Less massive objects are born from peaks with 
smaller height, v, and they acquire more ordered, h, and 
random, j, angular momenta. The density profile shape in 
the inner regions of haloes are particularly sensitive to an- 
gular momentum. In case of pure radial orbits, the core is 
dominated by particles from the outer shells. The larger is 
the angular momentum of particles, the longer they remains 
closer to the maximum radius. This give rise to a shallower 
density profile. Particles with smaller angular momentum 
will enter the core with a smaller radial velocity, with re- 
spect to particles in a pure radial SIM. Some particles have 
an angular momentum so large to be unable to fall into the 
core. So, particle with larger angular momentum will not 
reach the halo center and will not contribute to the cen- 
tral density, and this gives rise to flatter density profiles. 
Summarizing, larger ordered angular momentum (stronger 
interaction between the proto-dwarf and environment) gives 
rise to flatter density profiles, while a larger value of mass 
produces an increase of the central density contrast, and a 
steeper profile. 

Fig. Id shows the case of zero angular momentum and 
no baryons, for a halo of !O lo Af0. The solid line is the result 



of the model of the present paper, while the dashed line is the 
Einasto profile, characterized by a logarithmic slope varying 
continuously with radius (Navarro et al. 2004, 2008), and 
having a functional form: 

ln(p/p- 2 ) = (-2/ V )[(r/r- 2 y -1] (14) 

where the characteristic radius, r_2, is obtained by 
din p/dlogr — —2, p_2 = p(r-z), and 77 is a free param- 
eter, which for r\ ~ 0.17 reproduces fairly well the radial de- 
pendence of density profiles. In Fig. Id, we did not add the 
plot of NFW profile since recent simulations (e.g., Navarro 
et al. 2004, Navarro et al. 2010) shows that density pro- 
files from N-body simulations are better approximated by 
Einasto profiles. The correlation between steepness of den- 
sity profiles and angular momentum is very important to 
have a better understanding of LSB galaxies. We recall that 
LSB galaxies are more angular-momentum-dominated than 
normal galaxies of the same luminosity (McGaugh & de 
Blok 1998; Boisier at al. 2003), implying a flatter density 
profile with respect to giant galaxies. Another parameter 
which influence the final density profile is the baryon con- 
tent of the halo, namely the ratio of the baryon mass to 
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Figure 2. DM halos shape changes with baryon fraction. Same as previous figure, Fig. 1, but in panel (a) we reduced the value of 
baryon fraction of Fig. la (h,, fd,) to /d„/3, and in panel (b) we reduced the value of baryon fraction of Fig. lc to fd,/3; in panel (c) 
we increased the value of baryon fraction of Fig. la (h,/2, fd,) to 3fd,, and in panel (d) we increased the value of baryon fraction of 
Fig. lc to 3f d , . 



the total mass. In the following, we will study how changes 
in baryon content affects the density profiles. To this aim, 
in Fig. 2a, we plotted the same density profiles of Fig. la 
but with a baryon fraction fd, /3, and the value of specific 
angular momentum as left unchanged to h*. The effect of 
reducing the baryon content has the effect of steepening the 
density profiles, with respect to the reference case. The in- 
ner slope for haloes having the same mass as the reference 
case (10 s A/ Q , 10 9 M Q , 10 10 M Q ) are, respectively, given by 
q ~ 0.13, a ~ 0.16, and a ~ 0.32. Fig. 2b represents the 
haloes of Fig. lc (h,/2, fd,) but reducing the baryon frac- 
tion to fd, /3. The inner slope are steeper than the one seen 
in Fig. lc, namely we have a ~ 0.4, a ~ 0.6, and a ~ 0.7. 
Fig. 2c represents the same density profiles in Fig. la (ft*, 
fd, ) but in this case the baryon fraction is larger, namely 
3/d„ ■ Fig- 2d represents the same density profiles in Fig. lc 
(h,/2, fd,) but in this case the baryon fraction is larger, 
namely 3fd, ■ A larger amount of baryons imply a flattening 
of the density profiles. Fig. 2c shows density profiles with 
a ~ 0, for the three previous mass values, and Fig. 2d shows 
density profiles with a ~ 0.3, a ~ 0.45, and a ~ 0.51. The 
tendency to have flatter profiles with increasing baryon con- 
tent and shallower with larger baryon content is due to the 
fact that when more baryons are present the energy and an- 
gular momentum transfer from baryons to DM is larger and 
DM moves on larger orbits reducing the inner density. 



In Fig. 3, we check the model against G10 SPH simu- 
lations of dwarf galaxies in a ACDM framework. The figure 
was obtained following the prescription of G10 to calculate 
density profiles. In SPH simulations, G10 studied the for- 
mation of dwarf galaxies in a ACDM cosmology including 
baryonic processes, as gas cooling, heating from the cosmic 
UV field, star formation and SN driven gas heating was in- 
cluded. Dark matter particle masses in the high resolution 
regions had 1.6 x 10 4 Mq, and the force resolution, i.e., the 
gravitational softening, is 86 pc. They studied two dwarf 
galaxies, with mass 3.5 (DG1) and 2.0 x 10 10 M Q (DG2) in 
a ACDM cosmology with erg = 0.77, Q M = 0.24, Ov = 0.76, 
h = 0.73, and n b = 0.042. 

Fig. 3 plots Governato's DG1 (solid blue line) and DG2 
(solid black line) density profiles. Our profile is the dot 
dashed line. The plot shows a good agreement of our density 
profile with SPH results. 



3.2 Comparison with observed galaxies 

In the following, we use mass models for the dark matter 
component of different dwarf galaxies, from a dwarf with 
flat inner profile (NGC 2976), to an intermediate one (NGC 
5949) to one with a steep profile (NGC 5963) and compare 
these with those obtained from the model of Sect. 2. 
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3.2.1 NGC 2976 

In Figs. 4a-b, we apply the model of Sect. 2 to find density 
profiles, as done in sect. 3.1, that we then compare to the 
minimum disk and maximum disk rotation curve of NGC 
2976, respectively. Before, we summarize the characteristics 
of the quoted dwarf fully described in Simon et al. (2003) 
(hereafter S03). NGC 2976 is a regular Sc dwarf galaxy, a 
pure disk system, unbarred, and bulgeless (see Fig. 1 of S03). 
The dwarf, located in the group of M81 has a total mass, 
at 2.2 kpc, of 3.5 x 10 9 M o . S03 obtained, for the stellar 
disk, average values of 0.48 ± O.O2M0 /Lqk in K band and 
1.O7±O.O7M0/L0_r in R band. Stellar disk and gaseous disk 
constitute the fundamental reservoir of baryons of the dwarf. 
NGC 2976 has also a molecular disk and a gas disk, whose 
mass was estimated to be 1.5 x 10 s M©, by Appleton, Davies, 
& Stephenson (1981), and Stil & Israel 2002. The kinematics 
of the gaseous disk is highly complex, exhibiting large non- 
circular motions near the center which could be he reflection 
of a triaxial halo (Kazantzidis et al. 2010). S03 determined, 
for NGC 2976, the minimum and maximum disk rotation 
curve, represented by the black circles with error-bars in 
Fig. 4a, 4b, respectively. S03 also fitted the rotation curve 
with a power-law from the center to r ~ 110" (1.84 kpc) 
(solid line in Fig. 4a) and they found that the correspond- 
ing density profile, representing the density corresponding to 
the total mass (baryon plus dark matter), is well fitted by 
ptot = 1.6(r /lpc)-°- 27±0 - 09 Me/pe 3 . In Fig. 4a, the dashed 
line represents the result of the model of section 2 of the 
present paper considering the total mass (sum of dark mat- 
ter and baryons) (note that this is the only case in the paper 
where we plot he total mass: dark matter and baryons in or- 
der to compare results with the minimum disc rotation curve 
of S03), using A ~ 0.04, compatible with the value of A given 
in Munoz-Mateos et a. (2011). The baryon fraction in the 
case of NGC 2976 could be obtained as (M gas + M«)/M to t, 
where M tot = 3.7 x 1O 9 M (S03), M gas = 1.5 x 10 8 M Q (5% 



of the total mass) and M* = 5 x 10 8 M Q (14% of the total 
mass) (S03). However, the values of the mass given refers, as 
previously reported, only to the inner 2.2 kpc of the dwarf, 
so the quoted ratio would give just the baryon fraction of 
the inner part of the dwarf. For NGC 2976 as also for NGC 
5949, and NGC 5963 there is not enough data to determine 
the total mass and then the baryon fraction, since the dark 
matter halo is vastly more extended than any of the baryon 
tracers that are needed to measure the mass. So, one have 
to make a series of very broad assumptions, as done in Mc- 
Gaugh et al. (2010). In the case of NGC 2976, we know the 
inner masses and the circular velocity V c — 75 km/s. Gov- 
ernato et al. (2007) simulated a disk galaxy with V c = 70 
km/s and found that it has a virial mass of 1.6 x 10 n M s , 
similarly to Kaufmann et al. (2007), who simulated galaxies 
with V c - 74 km/s getting M vir = 1.1 x 10 n M Q ). So, as 
for Fig. 1, 2, we shall use McGaugh estimates, taking ac- 
count of the known rotation circular velocity of NGC 2976 
and the simulated mass in Governato et al. (2007). This 
gives a value fd — 0.1. The plot shows that, on the entire 
length of the data, the model of this paper gives a very 
good fit to the rotation curve, better than the power-law 
fit obtained by S03 or a NFW fit (here not plotted). The 
shape of the density profile of the dark matter halo (Fig. 
4b), was obtained by S03 removing the rotational velocities 
contributed by the baryon components of the galaxy, ob- 
taining the dark matter density profile of the maximal disk, 
Pdm = O.l(r/lpc)~ ool±ol3 M /pc 3 In Fig. 4b, the black 
circles with error-bars represent the dark matter rotation 
velocities for the maximal disk, as previously reported, and 
the solid line is the power-law fit (for 14" < r < 105") ob- 
tained by S03. Now, from the previous discussion, the slope 
of the total density profile of the galaxy (a dm ^ 0.27±0.09), 
composed of dark matter stars and gas, represents the ab- 
solute upper limit for the slope of the dark matter density 
profile. The lower limit to the dark matter density profile, 
obtained from the maximal disk, is aoM ^ 0.01 ± 0.13. 
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S03 concluded, after excluding strong and very young star- 
burst, that the dark matter density profile is bracketed by 
Pdm = r _0 - 17±0 - 09 and poM = r°. In Fig. 4b, the dashed 
line represents again the result of the model of section 2 of 
the present paper, now considering only dark matter. Also in 
this case, the model of this paper gives a very good fit to the 
rotation curve, better than NFW model and the power-law 
fit. 

In the previous discussion, we saw that the den- 
sity profile of NGC 2976 is very shallow, pnu = 
0.1(r/lpc)~ ' 01±ai3 M Q /pc 3 , compatible with some results 
in literature (e.g., de Blok et al. 2003), but in contradic- 
tion with the N-body simulations results (Stadel et al. 2009 
found a minimum value of the slope a = —0.8 at 120 pc). 
As discussed by several papers (e.g. DP09; Governato et 
al. 2010), the shallow central density profile of NGC 2976 
does not necessarily imply a problem for CDM model. Here, 
we only recall that one reason for the discrepancy between 
simulations and observations is due to the fact that simula- 
tions neglect the effects of the baryons on the dark matter 
halo. In the peculiar case of NGC 2976, baryons dominate 
the central region of the galaxy out to 220 pc, for the lower 
limit to mass-to-light ratio, and out to a radius of 550 pc, 
in the case of the maximal disk (S03). So, summarizing, the 
flat density profile of NGC 2976 is produced, as discussed in 
DP09, by the role of angular momentum, dynamical friction, 
and the interplay between DM and baryon component. The 
importance of tidal torques in shaping dwarf galaxies den- 
sity profile has been also suggested by several authors (e.g., 
Stoehr et al. 2002, and Hayashi et al. 2003). According to 
the previous works, smaller satellites of dark matter halos 
suffer tidal stripping, and tidal effects that mainly strip mass 
from the external regions of orbiting object and produce an 
expansion of the halo with the effect of a decrease in the 
central density at each peri-centric passage. So, even if the 
effect of tidal field is more evident in the outer part of the 
dwarf, also the inner profiles is influenced and it is shallower 
than their original profiles. However, tidal stirring is more 
efficient in forming dwarf spheroidal (dSph) galaxies start- 
ing from progenitors that are late type dwarfs affected by 
tidal forces (Kazantzidis et al. 2010, Lokas et al. 2010, Sales 
et al. 2010). The mechanism now quoted could flatten NGC 
2976 density profile if it is member of M81 and if it is at such 
a distance that tidal effects of M81 can influence it. If we 
assume for M81 the total mass determined by Karachentsev 
et al. (2002), namely 1O 12 M , in order that M81 tidal field 
is comparable to NGC 2976 gravity the two objects should 
approach each other at a distance of 20 kpc. Nowadays, their 
projected distance is known to be 79 kpc, and this means 
that, at present-time, NGC 2976 is probably unaffected by 
the interaction. This does not imply that if nowadays NGC 
2976 does not seem to be participating in the tidal inter- 
action currently taking place between M81, M82, and NGC 
3077 (Yun, Ho, & Lo 1994), that it was not interacting with 
them in the past. The existence of a HI streamer connecting 
M81 and NGC 2976 (Appleton et al. 1981) is probably the 
rest of a past interaction. More recently, Boyce et al. (2001) 
showed that the gas is composed of a single tidal bridge con- 
necting the two galaxies (see their Fig. 2a). In summary, af- 
ter the proto-structure decouples from the expanding back- 
ground and turns around, the growth of angular momentum 
is reduced to a second order effect, but subsequent interac- 



tions (e.g., the one between M31 and NGC 2976) have an 
important role in the overall evolution of the total angular 
momentum of the protostructure. This is confirmed by the 
results obtained from interacting systems, where the spin 
is visibly perturbed by the presence of a nearby compan- 
ion. Tidal interaction of a dwarf with a neighbor in later 
phase of evolution of the dwarfs can increase (but see also 
Cervantes-Sodi 2010a) the effects of tidal torquing in the 
phase of proto-dwarf which furnish the angular momentum 
to the dwarf itself. In the present paper, we did not consid- 
ered the tidal interaction among M31 and NGC 2976, which 
would further shallow its density profile. 

3.2.2 NGC 5963 and NGC 5949 

In Figs. 5a, b, we studied two other galaxies, NGC 5963 
and NGC 5949. NGC 5963 is a SC galaxy in the NGC 
5866 group (Fouque et al. 1992). At a projected distance of 
430 kpc, there is another galaxy of the group, namely NGC 
5907, distance so large to lead us think that NGC 5963 is 
not interacting with NGC 5907, and similar considerations 
lead to conclude that it is currently not interacting with 
its neighbors. NGC 5963 has a similar blue magnitude as 
NGC 2976 but has a mass four times larger (1.4 x 10 10 M Q ). 
From the center outside, the galaxy is characterized by a 
bar-like structure, outside a disk-like region and a LSB, ap- 
proximately exponential disk. Usually, it is assumed that the 
central region of NGC 5963 is a bulge at the center of a very 
faint disk (e.g., Kormendy & Kennicutt 2004). Non-circular 
motions are present in the galaxy (S03). 

NGC 5949 is a SBc galaxy with similar profiles of sur- 
face brightness to those of NGC 2976 (S03), even if NGC 
5949 has a larger blue magnitude with respect to NGC 2976, 
namely Mb = —18.2, and also a larger mass (M ~ 1O 1O M0). 
Another similitude between the two objects is that both have 
a nucleus, a shallow (steep) inner disk (outer disk). 

The procedure to obtain rotation curves and dark mat- 
ter density profiles of the two dwarfs is the same to that 
of NGC 2976 (S03). The limits placed by S05 on a DM are 
aoM = 0.79 ± 0.17, odai = 0.93 ± 0.04 for the maximum 
and minimum disk of NGC 5949 and a DM - 0.75 ± 0.10, 
ctDM = 1-41 ± 0.03 for the maximum and minimum disk of 
NGC 5963. The dark matter density profile slopes derived 
after subtracting the stellar disks are odm = 0.88 for NGC 
5949, and a DM = 1.20 for NGC 5963, respectively. In Fig. 
5a, the black points are the dark matter halo rotation curve, 
obtained by S05, after subtracting the stellar disk. The thick 
green, the dotted cyan, the short-dashed magenta, show, re- 
spectively, the power-law, NFW, and pseudo-isothermal fits 
obtained by S05, respectively. The long-dashed yellow curve 
is the dark matter halo rotation curve obtained from the 
model in the present paper. NGC 5963 is characterized by a 
very steep central density profile. S05 showed that a power- 
law fit with ct DM = 1.20 and a NFW fit with c = 14.9 and 
r s ~ 11, represents very well the rotation curve. A pseudo- 
isothermal model gives a much worse The long-dashed 
yellow curve obtained with the model of this paper is an 

5 Differently from Fig. Id, in Fig. 5a, b we compared the results 
to the NFW profile to be consistent with S05, who compared his 
rotation curves to the NFW profile. 
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Figure 4. Panel (a). Minimum disk rotation curve of NGC 2976. Black circles with error- bars represent the rotation velocities relative 
to the minimum disk obtained by S03. The solid line is a power-law fit to the rotation curve which corresponds to a density profile of 
p oc r -0,27 , obtained by S03. The dashed line plots the result of the model of the present paper, (b) Maximum disk rotation curve of 
NGC 2976. Similarly to panel (a), but in this case p <x r~ 01 . 



excellent fit. In order to obtain the quoted good fit, inside 
~ 100 pc we not also had to reduce the magnitude of h 
but also that of random angular momentum, j, (similarly 
to what done in Williams et al. 2004). In fact, as shown 
in Fig. Id, one obtain the steepest slope in the case the 
proto-structure is characterized by h — and baryon frac- 
tion Ft = 0, as in dissipationless N-body simulations. The 
density profile is approximated by an Einasto profile. This 
last has a slope of ~ —0.8 at ~ 100 pc and ~ —1.4 and 
at ~ lkpc. If the proto-structure contains baryons and the 
ordered angular momentum is not zero the density profile 
will be flatter (as seen in Fig. 1, and Fig. 2). In the case of 
NGC 5963, the baryon fraction is obtained as in the case 
of NGC 2976, by means of Mcgaugh et al. (2010), for an 
objet having V c — 120 km/s, and is given by fd — 0.12. 
As always reported, a good fit to NGC 5963, is obtained 
with a value of specific angular momentum h*/2 and reduc- 
ing the random angular momentum j to j/2 (see Fig. 6 of 
Williams et al. 2004), which steepens the density profile. In 
Fig. 5b, we plot the disk-subtracted rotation curve of NGC 
5949. Symbols are as in panel Fig. 5a. A power-law fit with 
cidm = 0.88, or a PI density profile, give good fits to the 
density profile, while a NFW model fit has the problem that 
fit parameters are not usefully constrained (S05). The long- 
dashed yellow curve obtained with the model of this paper 
is an excellent fit. In this case, V c — 80 km/s, and fd — 0.1 
and the value of specific ordered angular momentum /i*/1.6. 
Before going on, we want to recall that S03 and S05 analy- 
ses method differs in a few aspects from the other precedent 
cited studies. Firstly, their inner slope values are derived 
from a single power law fit to the entire rotation curve. Sec- 
ondly, most of their models do not take into account the 
gas component due to a lack of HI observations. However, 
the stellar disk has much more mass than gas, at the center 
of the studied galaxies. As shown by Bolatto et al. (2002), 



and S03, if one takes into account also the gas, one obtains a 
change in M/L of ~ 20%. The slope of the dark matter den- 
sity profile decreases of a factor 2%-12%, if one arbitrarily 
increases the stellar mass-to- light ratio (M*/L) to simulate 
this effect. As shown by S05 (section 2.2), if one allows that 
stellar disk have a non-zero thickness, one obtains an offset 
of this decrease. 

Finally, I would like to add that the same kind of pro- 
cesses operating in the late-type dwarfs that we previously 
discussed take place in LSBs, and that LSBs case was al- 
ready studied in DP09 (Fig. 4). 



3.2.3 Explaining the differences 

The quoted examples show that not all dwarfs are character- 
ized by flat inner cores (like NGC 2976), some have interme- 
diate slopes between pseudo-isothermal profiles and NFW 
profiles (NGC 5949) other are well fitted by NFW profile 
(NGC 5963). One question we could ask is: why NGC 5949, 
NGC 5963 (and also NGC 4605; NGC 6689) have steeper in- 
ner profiles, in some case (NGC 5963) compatible with the 
NFW model? In Sect. 3.1, we have seen how larger values 
of the baryons fraction and tidal interaction with neighbors 
produces flatter profiles. However, Fig. 1 and Fig. 2 show 
that an increase in angular momentum has a role more im- 
portant than an increase in baryon fraction (increase of a 
factor 2 in angular momentum produces similar effects than 
increase of a factor 3 in baryon fraction). For what concerns 
baryon fraction, the baryon fraction in the case of NGC 5963 
is fd ~ 0.12 and it is similar in the case of the other two 
galaxies (NGC 2976, NGC 5949). One should take into ac- 
count that it is not only important if the baryons fraction 
is similar for different dwarfs (as in the case of NGC 2976, 
NGC 5949, NGC 5963) to have similar density profiles (all 
the rest being the same), but is very important how it is 
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Figure 5. Panel (a). Disk-subtracted rotation curve of NGC 5963. Similarly to Fig. 4, the black points represent the dark matter 
halo rotation curve, obtained by S05, after subtracting the stellar disk. The thick green, the dotted cyan, the short-dashed magenta, 
show, respectively, the power-law, NF W, and pseudo-isothermal fits obtained by S05, respectively. The long-dashed yellow curve is the 
dark matter halo rotation curve obtained from the model in the present paper. Panel (b). Disk-subtracted rotation curve of NGC 5949. 
Symbols are as in panel (a). 



distributed in the galaxy (e.g., Hoeft et al. 2006, show in 
their figure 5 how the baryons fractions can strongly change 
going from the center to the outskirts of the dwarf). Now, 
the inner baryon fraction of NGC 2976 and NGC 5963 (us- 
ing S03, S05 data and Bosnia et al. 1988 data), at 2 kpc, are 
similar, and of the order of /;,. This very high inner baryon 
fraction are due to a funneling of baryons toward the center 
of the dwarfs. In fact, in the case of NGC 2976, Williams 
et al. (2010) results suggest that star formation is shutting 
down from the outside-in. It is possible that depletion of 
gas from the outer part of the disk, and formation of stars 
in the central part of the galaxy arose because of an inter- 
action of NGC 2976 with the core of M81 group more than 
1 Gyr ago. This last event, could have stripped halo gas and 
then produced an inflow of it from outer disk to the cen- 
ter of galaxy (Williams et al. 2010). So, in NGC 2976, the 
funneling of baryons toward the central region produces a 
larger exchange of angular momentum through baryons and 
dark matter and an inner flatter profile with respect to other 
dwarfs. We know that the surface brightness profile of NGC 



5949 is very similar to that of NGC 2976, and probably it 
has similar inner baryon fraction of NGC 2976. For what 
concerns, NGC 5963, differently from the other galaxies, in- 
ner and outer surface brightness profiles are characterized 
by a sharp transition. As a consequence, the derived stel- 
lar disk rotation curve have a peak at small radii followed 
by a steep fall, with the consequence that the galaxy has 
outer parts dominated by dark matter (see Fig. 6b of S05). 
As previously discussed, the galaxy could contain a pseudo- 
bulge, and so like in NGC 2976 some process have funneled 
baryons toward the center of the galaxy. This means that in 
the case of these three dwarfs the baryon fraction is playing 
a similar role in shaping the density profile. According to our 
model, another factor in shaping density profile is the mass 
of the dwarf. Among the galaxies studied by S05, namely 
NGC 2976, NGC 4605, NGC 5949, NGC 5963, and NGC 
6689, NGC 5963 has the highest central surface brightness 
and the most unusual surface brightness profile. It does also 
have the highest mass, in the sample and the highest mag- 
nitude (Ali = —19.1). From the point of view of our model, 
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a larger mass imply a steeper profile (however in the case of 
NGC 5963, NGC 5949, and NGC 2976 this difference is very 
small (see Table 7 of S05). Also in a recent paper, dB08, by 
using the THINGS sample, showed that in smaller objects 
core dominated halo is clearly preferred over a cusp-like halo, 
while for massive, disk dominated galaxies, PI and NFW fit 
apparently equally well. In other terms, for low mass galax- 
ies a core dominated halo is clearly preferred over a cusp- 
like halo, while for massive, disk dominated galaxies, all halo 
models give fits of similar quality. This because more massive 
galaxies tend to have a more compact light-distribution, and 
therefore a more compact stellar mass distribution as well. 
This makes the contribution of the disk to the inner rotation 
curve critically dependent on the stellar mass-to- light ratio, 
and the initial mass function (IMF). So one could imagine a 
situation where the M/L* could be changed to produce flat 
profiles for massive galaxies, and these changes do not have 
to be very large. Low mass galaxies have instead a reduced 
sensitivity to the M/L* values, implying that changes of the 
quoted ratio does not influence the profile shape. 

For what concerns, the role of tidal torques, we know 
that in the case of NGC 2976 we need a higher value of 
the specific angular momentum [h, A = 0.04) than in the 
case of NGC 5963 and NGC 5949 to fit the profile. Tidal 
torques influence not only dwarfs during their formation, but 
if the environment contains close galaxies, tidal torques can 
influence the shape of the density profile. A search in NED 
(Nasa Extragalactic Database), concerning NCG5963 and 
the other quoted galaxies, gives to us the following results: 
as previously reported, NGC 2976 is a satellite of M81, NGC 
5949 closest neighbor is a dwarf galaxy 300 kpc away. NGC 
5963 is part of a small group; the closest neighbor is a dwarf 
galaxy 200 kpc away, and the massive galaxy NGC 5907 
(3x 10 10 M o , Lequeux et al. 1998) is about 400 kpc away. So, 
NGC 5963 is the only other galaxy that could be regarded 
as having a similar environment to NGC 2976, even if the 
actual projected distance between M81 and NGC 2976 is 
79 kpc and in the case of NGC 5963 the neighbors are at 
a much larger distance. Moreover, although for NGC 2976 
there is good evidence that it has interacted with M81, there 
is not really significant evidence for an interaction between 
NGC 5963 and its neighbors. In other terms the galaxy can 
be considered isolated and, as seen in Sect. 3.2, this imply 
that one cannot expect flattening of the density profile by 
tidal interactions, like for NGC 2976. Similar conclusions are 
valid for NGC 5949. Differently, in the case of NGC 2976, the 
interaction with M81 further contributed to the flattening 
of the density profile. 

In an attempt to have a deeper insight on the issue (in- 
ner profile shape and environment), we have analyzed the 
rotation curves of a sample of dwarfs, LSB, and three galax- 
ies from THINGS. We fitted the rotation curves of the dark 
halo with power laws, as did by S05 and OhlO and related 
each galaxy to the so-called tidal or galactic isolation index, 
TI, whose negative values correspond to isolated^ galaxies 
of the general field and positive values to members of the 
groups. The quoted index was introduced by Karachentsev 



6 Isolated galaxies are galaxies located in low density environ- 
ments and selected according to fixed criteria (see Karachentscva 
1973; Karachentsev & Petit 1990; Vavilova et al. 2009) 



& Kashibadze (2006) in their study of the masses of the Lo- 
cal Group and of the M81 group, obtained by distortions in 
the local velocity field. It is defined as: 

8 = max[log(M k /Dl)} + C (15) 

where (Mk',Dk) are the mass and distance to the chosen 
galaxy of the k - th system and C is a normalization con- 
stant. Estimating O for a given galaxy is actually quite dif- 
ficult and typically possible only for systems in the Local 
Group. A cross match between our sample and the Cata- 
log of Neighboring Galaxies (Karachentsev et al. 2004) and 
Karachentsev and Kashibadze (2006) allows us to infer O 
for the sample. 

The result is plotted in Table 1. The first column iden- 
tifies the galaxy, the second the type of galaxy (dwarfs 
(type=0), LSB (type=l) and the three larger galaxies from 
THINGS (type=2), and the third the paper from which it 
was taken, last column is the so called tidal index, ©, defined 
as (Karachentsev et al. 2004; Karachentsev and Kashibadze 
2006). 

In the case of NGC 2976, we have TI — 2.7, indicating 
that the galaxy is non-isolated, and as previously seen its 
inner slope has a ~ 0.01. In the case of NGC 4605, we have 
TI = —1.1, indicating that the galaxy is isolated, and the 
inner slope is steeper (a ~ 0.78). In the case of NGC 2366, 
we have TI = 1, indicating that the galaxy is non-isolated, 
and the inner slope is ot ~ 0.32 (OhlO). As reported in 
the introduction, one of the most recent studies on mass 
modeling for dark matter component of galaxies is the work 
of OhlO, in which all the seven dwarfs chosen from THINGS 
are better described by core-like models. These galaxies are 
included in Table 1. However, we want to point out that at 
least 5 of the quoted dwarfs are characterized by a positive 
TI, namely NGC 2366 (T = 1), ICG2574 (T = 0.9), HoII 
(T = 0.6), DD053 (T = 0.7), HOI (T = 1.5), indicating 
that the objects are non-isolated, and one should expect 
flat profiles, as observed. It would be interesting to extend 
the previous set of OhlO to include isolated objects and to 
compare the density profile of the last with those of non- 
isolated dwarfs. 

From Table 1, is evident that, in the limit of the small 
sample considered, a correlates with O. Moreover, the large 
value of the Spearman rank correlation coefficient, being 
C(a;0) = —0.81 (for dwarfs case), strongly suggests that 
the trend observed is not an artifact due to small statistics. 
We want to add that we are working to build up a big- 
ger sample (also using the WHISP sample by Swaters et al. 
2002), and and make fits with power laws, and generalized 
NFW models. In order to further test the possible correla- 
tion between the inner slope and the environment, we will 
resort to two other commonly used indicators. The first, de- 
noted as Eio, is the projected number density of galaxies 
inside the circle of radius equal to the projected distance of 
the tenth nearest neighbor (Dressier et al. 1980), while the 
second, denoted as -D750, is the projected number density of 
galaxies within a circle of radius equal to 750 kpc centred 
on the system of interest (Treu et el. 2009). £10 and -D750 
positively correlate with O so that we can use them as sup- 
plementary tools to check whether the environment impacts 
the inner slope of the density profile. 

The quoted example and the discussion of the previous 
sections, indicates that tidal interaction, baryon fraction, 
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Table 1. Dwarf (type = 0), LSB (type = 1), and normal spirals 
(type=2) sample and constraints on model parameters. 



Id 



Type 



Ref 



(-) 



Q2+ - 08 
n SQ+ 06 

2 + 0.06 
-0.06 

3+0.06 
-0.06 
3 +0.07 
-0.07 
,+0.08 
-0.08 
-+0.07 
-0.07 
,+0.06 
-0.06 
,+0.06 
-0.06 
2+0.06 
-0.06 
7+0.05 
-0.05 

3 + 0.06 
-0.06 

2+0.06 
-0.06 
3+0.07 
-0.07 
3 + 0.1 
-0.1 
2+0.06 
-0.06 

-0.1 



0.78Z 
0.39j 
0.79j 
0.8lj 
0.75j 
0.51^ 
0.46^ 
0.38j 
0.07j 
0.43^ 
0.38^ 
0.09^ 
0.32^ 
0.38j 
0.01" 1 



UGC1281 1 dBB02 -1.2 

UGC8490 Swll -1.1 

NGC4605 S05 -1.1 

M81dwB OhlO -0.9 

NGC5023 1 dBB02 -0.5 

NGC4736 2 dBB02 -0.5 

NGC3274 1 dBB02 -0.3 

NGC2403 2 dBB02 

UGC7524 Swll 0.1 

UGC7559 Swll 0.1 

NGC7793 2 dBB02 0.1 

HoII OhlO 0.6 

DD053 OhlO 0.7 

IC2574 dBB02 0.9 

NGC2366 dBB02 1.0 

Hoi OhlO 1.5 

NGC2976 dBB02 2.7 

Swll (Swaters et al. 2011) 



and environmental effects, influence the slope of the dwarfs 
density profile. This does not come unexpected, since we 
know that on the other hand, the morphology-density re- 
lation illustrates that environment is an important factor. 
The existence of a correlation between galaxy morphology 
and local density environment (Dressier 1980, Goto et al. 
2003, Park et al. 2007) is a motivation to search for cor- 
relations between local environment and internal properties 
of galaxies. One of this correlation of interest to us, is the 
possible correlation between A and environment. Concern- 
ing this issue, there are different opposite results. Maccio 
et al. (2007) found an absence of correlation between A and 
environment in cosmological simulations, although the oppo- 
site was found by other studies using similar techniques (e.g. 
Avila-Reese et al. 2005, Cervantes-Sodi et al. 2008a,b, 2010a, 
2010b). Avila-Reese et al. (2005) employed a ACDM N-body 
simulation to study the properties of galaxy-size dark mat- 
ter haloes as a function of global environment, where A was 
one of the studied properties of the halo. Cervantes-Sodi 
found a marked correlation of A with mass (Cervantes-sodi 
et al. 2008a, b), a correlation between A and cluster en- 
vironment (Cervantes-Sodi et al. 2010b), and showed that 
events such as interaction with close neighbors play an im- 
portant role in the value of the spin for the final configura- 
tion (Cervantes-Sodi et al. 2010b). Even the spin alignment 
of dark matter haloes in different environments, has been 
studied in some N-body simulations (Aragon-Calvo et al. 
2007). From the observational point of view, there has not 
been done much study on the effects of environment on the 
slopes of the dark matter halos, simply because there are 
not many large samples needed for that, yet. As previously 
reported, there is a first sample (WHISP by Swaters et al. 
2002) and to this will be added in the next few years, the 
LITTLE THINGS (PI: D. Hunter) and VLA-ANGST (PI: 
J. Ott), and they should give access to high-resolution HI 



observations from ~ 50 to 60 dwarf galaxies in various en- 
vironments. Moreover, THINGS in that respect is a sample 
that was specifically selected not to contain too many inter- 
acting galaxies (apart from one or two show case galaxies), 
and is therefore not an ideal sample to look at the effects 
mentioned. Moreover, the techniques actually used to infer 
the rotation curve makes difficult to address the problem of 
the quoted correlation. The problem is that when a galaxy 
is undergoing an interaction, the orbits of the gas particles 
will change, especially those in the outskirts. So the entire 
assumption of the gas going around in circular orbits, and 
the analysis using the ROTCUR approach, will break down. 
One could, of course, still try to analyze the velocity field 
and, making some reasonable assumptions, come up with a 
rotation curve. But if one then find that the inner slope is 
smaller, there is the problem of how confident can one be in 
that result: the low value of the slope could be due to the 
fact that we are fitting a model of circular velocities to what 
in reality are non-circular motions. One partial solution to 
the problem would be that of studying dwarfs non strongly 
interacting with neighbors and to compare with dwarfs in 
voids or to wait for more sophisticated techniques. Clearly, 
many work is needed from the theoretical and observational 
point of view in order to have deeper insights on the role of 
environment and density profiles of dwarfs and other impor- 
tant questions require further research, including the mech- 
anisms governing dwarf galaxy evolution, and the origin of 
the morphology-density relation. 



4 CONCLUSIONS 

Using the SIM model introduced in DP09, we studied how 
the shape of density profiles of dwarfs in the mass range 
10 8 — 1O 1O M0 changes when angular momentum originated 
by tidal torques, and baryon fraction are changed. As a first 
step, we calculated the reference density profiles following 
DP09, plotted in Fig. la, then we calculated density pro- 
files of haloes having ordered angular momentum in a range 
of values compatible with dwarf galaxies haloes and baryon 
fraction as obtained by McGaugh et al. (2010). We found 
that density profiles steepen with increasing mass of the halo 
while increasing angular momentum the density profiles flat- 
ten, in agreement with previous works (Sikivie et al. 1997; 
Avila-Reese et al. 1998; Nusser 2001; Hiotelis 2002; Le Del- 
liou & Henriksen 2003; Ascasibar et al. 2004; Williams et al. 
2004) . In the case tidal torquing is shut down (zero angular 
momentum) and baryons are not present, the density profile 
is very well approximated by an Einasto profile. Increasing 
baryon fraction one gets flatter profiles. This is due to the 
fact that when more baryons are present the energy and an- 
gular momentum transfer from baryons to DM is larger and 
DM moves on larger orbits reducing the inner density. We 
then repeated the same calculations for three dwarfs stud- 
ied by S03, S05, namely NGC 2976, NGC 5949, and NGC 
5963. NGC 2976 has a core profile, NGC 5963 a cuspy pro- 
file and NGC 5949 an intermediate one. We calculated the 
baryon fraction using some assumptions and McGaugh et al. 
(2010) results and then fitted the rotation curves obtained 
by S03, S05 with the theoretical curves obtained from our 
model. NGC 2976 is characterized by a larger value of angu- 
lar momentum (A ~ 0.04) with respect the other two dwarfs 
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(A ~ 0.025, for NGC 5949; A ~ 0.02, and random angular 
momentum, j, 1/2 of those of the reference haloes (Fig. la), 
for NGC 5963) . Moreover, there are evidences of past tidal 
interaction between NGC 2976 and M31. Appleton et al. 
(1981) discovered a faint HI streamer stretching from M81 
to NGC 2976, implying an interaction between the two ob- 
jects, and observations of Williams et al. (2010) showed a 
process of outer disk gas depletion and inner disk star for- 
mation likely triggered by an interaction with the core of 
the M81 group more than 1 Gyr ago. This interaction could 
be another reason why the density profile of NGC 2976 is 
so flat. In the case of NGC 5949 and NGC 5963 there is 
no evidence of this kind of tidal interaction and the clos- 
est neighbors are at several hundreds kpc. The small sample 
that we used, shows a correlation between the a and the 
tidal index O. 

Concluding, dwarfs galaxies which acquired larger angu- 
lar momentum by tidal interaction in their formation phase 
and/or interacted tidally in later phases with neighbors, 
show a flatter density profile than those who had less tidal 
interactions. 
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